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Abstract. We obtain approximate analytic expressions for the critical value of the 
total angular momentum of a non-relativistic test particle moving in the Kerr geometry, 
such that it will be captured by the black hole. The expressions apply to arbitrary 
orbital inclinations, and are accurate over the entire range of angular momentum for the 
Kerr black hole. The expressions can be easily implemented in N-body simulations of 
the evolution of star clusters around massive galactic black holes, where such captures 
play an important role. 
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1. Introduction and summary 

The behavior of large numbers of stars orbiting around a central massive rotating black 
hole is a subject of great current interest (for recent reviews, see [HE])- Recently, Merritt 
et al. [3] used N-body simulations to construct and evolve a cluster of stars around a 
massive black hole such as the one at the center of our own galaxy, and studied the 
perturbations induced by those stars on the orbits of hypothetical stars that might be 
used to test the black hole no-hair theorems [HE]. In separate work, Merritt et al. 
also analyzed the rate of injection of stars onto highly eccentric orbits around a massive 
black hole [0]. Such orbits would lead to an "extreme mass-ratio inspiral" (EMRI) that 
would result in the emission of detectable gravitational radiation. 

Generally, such N-body simulations must be evolved for long periods of time in 
order to allow for all correlations and relaxation mechanisms to be established, so as to 
develop realistic estimates for rates of particular events, such as EMRI production. In 
the case of [6], the timescale was 10 Myr. 

The simulations of Merritt et al. [3] included the Newtonian gravitational attraction 
between the stars and the black hole and between each other, as well as post-Newtonian 
effects associated with the massive, spinning central black hole, including the dragging 
of inertial frames and the quadrupole moment of the black hole's field. The simulations 
of [6] restricted attention to a non rotating black hole, but included the capture of stars 
by the black hole. A star-star encounter can place one star on an orbit with sufficiently 
small angular momentum that it can be captured directly by the black hole. This is a 
non-trivial effect: in typical simulations carried out in |6j, out of 50 stars in the cluster, 
17 were captured by t = 2 Myr, and 30 were captured by t = 10 Myr. 

In [6] the condition for a capture was simple: if the total angular momentum per 
unit mass L of a star after an encounter with another star was smaller than the critical 
value L c = 4GM/c, where M is the mass of the black hole, the star was deemed to be 
captured, the mass of the black hole was augmented by the star's mass, and the hole 
was given a suitable recoil. This value of L c corresponds to the well-known angular 
momentum of an unstable circular orbit in the Schwarzschild geometry with conserved 
relativistic energy per unit mass E = 1. Grossman et a/.[7j call this the innermost bound 
spherical orbit (IBSO). The latter condition is appropriate here, because the stars are 
very far from the black hole in non-relativistic orbits when they encounter one another. 
Accordingly, one can treat them as being barely bound (E = 1) from the point of view 
of the black hole's strong gravity. 

However, there is every reason to expect that the typical massive black hole at 
the center of a galaxy, including our own, will be rotating, perhaps even very rapidly 
rotating. As a consequence, the capture of a star will not be isotropic: stars in prograde 
orbits will be less likely to be captured than stars in retrograde orbits. For example, for 
E = 1 stars with orbits restricted to the equatorial plane, the critical angular momentum 
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per unit mass for capture is given by 

L c = ^(2 + 2v / T^), (1) 

where a is the dimensionless Kerr parameter, related to the black hole's angular 
momentum J by a = Jc/GM 2 (0 < a < 1). The upper(lower) sign corresponds to 
prograde(retrograde) orbits (see, for example, [EJ |9]). This anisotropic capture could 
have important consequences for the evolution of the cluster, imparting a net angular 
momentum to it, for instance, because of the preferential loss of stars moving in a 
retrograde manner relative to the black hole. 

For orbits out of the equatorial plane, the problem is much more complicated, 
despite the presence of the additional "Carter" constant of motion, which makes the 
mathematical problem completely integrable. Generally one must resort to numerical 
solutions to study the behavior of orbits near this critical point [TOj [TTJ. 

Instead, we have found an analytic expression for the critical angular momentum 
for arbitrary orbits that is approximate, but surprisingly accurate over almost the whole 
range of black hole spins. It has the form 



GM I / 1,2,^ 

L c = I 2 + 2 y 1 — a cos % — -a 2 sin iF(a, cos i) I . (2) 

Because the total angular momentum is no longer a conserved quantity, we define L 
in terms of the Carter constant: L = \[C (all quantities are suitably scaled so as to 
be independent of the mass of the star). We also define cosz = L z /L, where L z is 
the conserved z-component of angular momentum per unit mass. In the Schwarzschild 
limit and the Newtonian limit, L is the total conserved angular momentum and % is 
the orbital inclination, with < % < tt/2 corresponding to prograde orbits (L z > 0) 
and tt/2 < i < n corresponding to retrograde orbits [L z < 0). The function 



F(a, cosz) = 1 + (a/2)cosz + ... is a series expansion in a shown below in Eq. (15). 
With the series carried to 0(a 4 ), the solution ^ agrees with numerical solutions for the 
critical angular momentum to better than 0.5 percent for < a < 0.9; and to better 
than 5 percent for 0.9 < a < 0.99. For the special case of the extreme Kerr black hole 
(5 = 1), we obtain a separate, but very accurate analytic fit for L c as a function of cosi. 
Equation ^ can be easily implemented in an N-body code as a capture criterion. 

The rest of this paper provides details. In Sec. [2] we review the basic equations for 
motion in the Kerr geometry, and in Sec. [3] we obtain the critical value of L for capture 
in the case of E = 1. Concluding remarks are made in Sec. |4j Henceforth, we use units 
in which G — c — 1. 

2. Basic equations for geodesic motion in the Kerr geometry 



The Kerr metric is given in Boyer-Lindquist coordinates by 

ds 2 = -II — dt 2 + ^rdr 2 + p 2 d6 2 — sin 2 Odtdcj) 

V p 2 / A p 2 
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r 2 + a 2 + sin 2 #d0 2 , (3) 

/r / 

where M is the mass and a is the Kerr parameter, related to the angular momentum 
J by a = J/M, with |o| < M required for the metric to describe a black hole; 
p 2 = r 2 + a 2 cos 2 8, and A = r 2 + a 2 — 2Mr. We will assume throughout that a is 
positive. 

Timelike geodesies in this geometry admit four conserved quantities: energy per 
unit mass of the particle E, angular momentum per unit mass L z , Carter constant C, 
and the norm of the four velocity u a = dx a /dr, given by 

E = -u = -g 00 u° - g H vfi , (4 a) 

L z = = g 0<f) u + g H u* , (46) 

C = p\u e f + sin" 2 6L 2 Z + a 2 cos 2 6(1- E 2 ) , (4c) 

gapvTu* = - 1 . (Ad) 

The version of the Carter constant used here has the property that, in the Schwarzschild 
limit (a — > 0), C — > L 2 , where L is the total conserved angular momentum. 
Solving Eqs. (4a) and (Ab) for u° and v?, and Eq. (4c) for u e , we obtain 

_ E9u + L z go4, 



u ' 



u 



Asin 2 ^ 

Ego<f, + L z g 00 



A sin" 9 ' 

p- 2 [C - sin" 2 9L 2 - a 2 cos 2 6(1 - E 2 )} 1 ' 2 . (5) 



Substituting these into Eq. (Ad) and solving for u r yields, after some simplification, 



\ r pA )\ (0 



where 



y r / \ A a 2 2Ma 2 \ n2 A/ C\ a 2 L 2 AMaEL z 

This potential can then be used to study turning points of the radial motion, and will 
be the key ingredient in studying the capture of particles. 



3. Critical angular momentum for capture 

We now want to find the critical values of the constants such that an orbit of a given 
energy E, angular momentum L z , and Carter constant C will not be "reflected" back 
to large distances, but instead will continue immediately to smaller values of r and be 
captured by the black hole. The turning points of the orbit are given by the values 
of r where u r = 0, thus where V(r) = 0. The critical values of E, C and L z are 
those for which the potential has an extremum at that same point, that is where 
d[(r 4 / p 4 )V (r)} I dr = 0. But since we are choosing a point where V(r) = 0, this is 
equivalent to dV(r)/dr = 0. The chosen sign for V(r) also dictates that this point 
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should be a minimum of V(r), that is that d 2 V(r)/dr 2 > 0, corresponding to an unstable 
extremum. 

Stars will be injected into orbits that could pass very close to the black hole via 
interactions with other stars far away from the black hole, in the non-relativistic regime. 
After the encounter, the orbit of such a star can be characterized by its Newtonian 
semi-major axis a (not to be confused with the Kerr parameter a), which will be large 
compared to M, and thus the energy per unit mass will be given approximately by 
E « (1 — M/2a) m 1. If this encounter can be regarded as the final significant 
perturbation of the orbit by other stars, this energy, along with the other constants 
of motion in the Kerr geometry (L z , C) will be conserved for the rest of the orbit. As 
a consequence, we shall focus attention on orbits for which E — 1. 

In this case, linear combinations of the conditions V(r) = and dV(r)/dr = 
yield two quadratic equations for r. These correspond to special cases of, for example, 
Eqs. (2.12) and (2.13) of [12], or Eqs. (Al) and (A2) of [7J. Solving for r from one 
equation and substituting back into the other equation yields an algebraic equation for 
the critical value of C in terms of L z , M and a. 

In order to make the transition to spherical symmetry and the connection to the 
Newtonian portions of the orbits more natural, we will define the new constants of 
motion, L > and cosi, related to C and L z by 



In both the spherical and Newtonian limits, L and i correspond to the total angular 
momentum and the orbital inclination, respectively. In the strong-field regime close to 
the black hole, these quantities remain constants of the motion, but their interpretation 
is less transparent. 

For orbits in the equatorial plane, L = \L Z \, and % = for prograde orbits or i = it 
for retrograde orbits. For general orbits, % ranges from to tt/2 for prograde orbits, and 
from 7r/2 to ir for retrograde orbits. 

Scaling out the factors of M by defining the dimensionless variables L = L/M, 
a = a/M, we find that the algebraic equation for the critical value L c takes the form 



where Q = \ j\L\(L\ — YlL 2 c + 24aL c cosi — 12a 2 ) 2 ]. The critical radius corresponding 
to this unstable circular orbit is given by 




(8) 



= Q{{\ - a 2 sin 2 i)L 8 c - Ad cos iL 7 c - 2 [8 - S 2 (3 + 7sin 2 i)]L; 
+ Ad cos i[2A - d 2 (l + 9 sin 2 i)]L 5 c 
- a 2 [240 - 192 sin 2 i - d 2 (l + 18 sin 2 i - 27 sin 4 i)}L 4 c 
+ 64S 3 cosi(5 - 2 sin 2 i)L 3 c - 48S 4 (5 - 4 sin 2 i)L 2 c 
+ 96a 5 cosiZ c - 16a 6 }, 
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c 



(9) 



a cos i) 2 — 8a 2 sin 2 i 



Li - 12{L, 



a cos i) 2 — 12a 2 sin 2 i 



(10) 



c 
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Figure 1. Critical angular momenta L/M for capture vs. Kerr parameter a/M, for 
values of cos? ranging from —1 to 1. 



Considering first the case of an equatorial orbit, with i = or i = it, the condition 
becomes 



(L 2 t 4L C + 4a) (L\ ± 4L C - 45) (L c =)= a; 



o = '.— : 1 --t c i — 1 ; . (ii) 

L2(L4 - 12L 2 C ± 24fiL c - 12fi 2 ) 2 

The only solutions that correspond to values of r outside the event horizon and to 
maxima of the potential V(r) are 

L c+ = 2 + 2\/l — a , i = (prograde) , 

L c _ = 2 + 2Vl + 5 , i = vr (retrograde) . (12) 

The corresponding critical radius is given by 

f c ± = (1 + VlTa) 2 = L c ±Ta. (13) 

For 5 = 0, the solutions merge to L c = 4 and r c = 4, corresponding to the well-known 
critical angular momentum in the Schwarzschild metric for the unstable circular orbit 
with E = 1. 

For arbitrary values of i, Eq. (|9| must be solved numerically. The numerical 
solutions are plotted in Fig. 1, showing the critical value L c /M vs. a/M for inclination 
angles ranging from to it. 

Alternatively, we shall seek an analytic, albeit approximate solution. To this end, 
we consider a solution for L r of the form 



L c = 2 + 2-1/1 — a cos i a 2 sin iF(a, cos i) , (14) 
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which has the correct behavior both when a = and when i = or i = ir for arbitrary 
a. We then solve Eq. ^ in a power series in a, obtaining 

F(a, cosi) = 1 H — a cos 2 H a 2 (7 + 13 cos 2 i) H a 3 cosi(23 + 5 cos 2 i) 

v ' ; 2 64 v 7 128 V ; 

+ t^S 4 (55 + 340 cos 2 i - 59 cos 4 z) + 0(a 5 ) . (15) 
2048 

This series solution agrees with numerical solutions of Eq. (JoJ) to better than 0.5 percent 
for < a < 0.9; and to better than 5 percent for 0.9 < a < 0.99. For 5 = 1, the critical 
value can be fit to better than 0.5 percent by the function 

f 1.575 + 2.276V1 - 1.044 cosi : -1 < cosi < 

L c = I (16) 

[2/cosz : a/2/3 < cosz < 1 . 

The second solution is actually an exact solution for the range indicated. 

Thus for a star in an orbit with E = 1 and a given cosz = L z /L, plunge will occur 
for L < ML C . When the star is far from the black hole, the geometry is at the same 
time approximately Schwarzschild and approximately Newtonian. In this case L is the 
total angular momentum; it can be expressed in terms of the Newtonian semimajor axis 
a and eccentricity e by L 2 = Ma(l — e 2 ) = Mr p (l + e), where r p = a(l — e) is the 
apparent pericenter distance inferred from the Newtonian orbit. For large eccentricity 
(e ~ 1), then, capture will occur when 

r P < r pc = ^M(L C± ) 2 . (17) 

For Schwarzschild, the capture pericenter radius is r c = 8M; for extreme Kerr on 
the equatorial plane, the capture pericenter radius is r c = 2M(11.6M) for prograde 
(retrograde) orbits. 

4. Concluding remarks 

Implementation of this capture condition in an iV-body code will depend on where 
along a given orbit the star's orbit elements are determined in order to find r p , i and L z . 
Presumably it is after the last important stellar encounter. If the star is still relatively 
far from the black hole, then one can argue that relativistic corrections will be small 



compared to the variations shown in Eq. (14). If the star is close to the black hole 



when this evaluation is made, then it may be necessary to be more careful in applying 



Eq. (14), for example by incorporating post-Newtonian corrections in the evaluation of 
the constants of motion for the orbit. Given values of position and velocity of the star 
following the encounter, this can be done in a straightforward manner. 

It would be interesting to try to generalize the capture condition to value of E less 
than unity, so as to handle more relativistic clusters. 
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